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We study a clean superconductor in the Hall configuration, in the framework of a purely dissipative 
time-dependent Ginzburg-Landau model. We find situations in which the order parameter differs 
significantly from zero in a set of islands that appear to form a periodic structure. When the pattern 
of islands becomes irregular, it moves in or against the direction of the current and a Hall voltage 
is found. Tiny differences in the initial state may reverse the sign of the Hall voltage. When the 
average Hall voltage vanishes, the local Hall voltage does not necessarily vanish. We examine the 
influence that several boundary conditions at the electrodes have on these effects. 
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I. INTRODUCTION 

The Hall voltage in superconductors exhibits a rich va- 
riety of behaviors.^ In the Meissner state there is no Hall 
voltage,^ but in the mixed state Hall voltage is present 
due to vortex drag.'^''' In some cases, the sign of the Hall 
voltage is opposed to what would naively be expected.^ 

We shall consider a thin rectangular superconducting 
sample. Let a magnetic field be applied in the z-direction 
and let a total current / flow in the a;-direction. The 
sample will be assumed to be sufficiently long in the z- 
direction, so that physical quantities will be independent 
of z. The current will be assumed to flow in the entire 
range — oo < x < oo, but only the segment Q < x < L will 
be superconducting. We denote the thickness of the sam- 
ple by d; the regions y < and y > d are taken as insulat- 
ing. We will study the current dependence of measurable 
quantities within the framework of the time-dependent 
Ginzburg-Landau model (TDGL).^ Among the diversity 
of formulations of TDGL, we consider the simplest: 

dti^ = ~- [(-zV - A)2 ^ + (1 - T) {\^\^ f , 

(1) 

dtA = (1 - T) Re [V5 (-iV - A)^]- k^V x V x A . (2) 

Here ip is the order parameter, t is the time, A is the elec- 
tromagnetic vector potential, 77 is the ratio between the 
relaxation times of ■0 and A, k is the Ginzburg-Landau 
parameter, T is the temperature and / a random "force" 
that simulates thermal fluctuations. The units are cus- 
tomary, as e.g. in Rcfs. 7-10. The gauge is chosen such 
that the scalar potential is zero. 

We shall see that this configuration leads to the ap- 
pearance of a phenomenon which, to my knowledge, has 
not been previously encountered: spots where {ipl is sig- 
nificant, whereas most of the sample is practically in the 
normal state. We call these spots "superconducting is- 
lands." The expression "superconducting islands" is usu- 
ally intended for regions of superconducting material sep- 
arated by thin insulating barriers, but in our case these 
islands will form spontaneously in a uniform material. In 



a loose sense, superconducting islands may be regarded 
as the opposite of vortices, but their length scale is much 
larger. In appropriate situations, these islands form a 
periodic pattern. We shall see that the rearrangement of 
these islands is related to the appearance of Hall voltage. 



II. SELECTION OF THE PROBLEM 

The situation we consider is as follows: the ap- 
plied magnetic field is kept fixed at 0.5ffc2(0) = 
0.25<I>o/7r^^(0), where $0 is the quantum of flux and ^(T) 
is the coherence length at temperature T. At this field, 
the sample is in the mixed state. Initially, there is no net 
current and the situation is static. Then the current in 
the x-direction is gradually increased, until the film be- 
comes normal. If the current increases sufficiently slowly, 
we may argue that we have a quasistationary situation 
and thus evaluate the properties of the superconductor 
as functions of the current. 

The boundary conditions that are usually assumed in 
the TDGL treatment are continuity of the magnetic field 
and {> ■ (— — A) ip — 0, where v is a vector perpendic- 
ular to the superconductor-insulator interface. However, 
this condition implies that the electric field is parallel 
to the interface, and I therefore suspected that it might 
not be appropriate for the study of the Hall voltage. In- 
stead, I used the refined boundary condition suggested 
in Eq. (8.26) of Ref. 9 at the boundaries y = and 
y — d. It turned out, however, that the same results 
are obtained without this refinement. At the electrodes 
(x — and x — L) I have considered two different bound- 
ary conditions. One case was that of periodic boundary 
conditions, which are frequently used to mimic an infi- 
nite sample in the x-direction; the other case was the 
Dirichlet condition, as appropriate for normal electrodes 
in which superconductivity is strongly suppressed. Peri- 
odic boundary conditions arc not physically justified, but 
they lead to results that are simpler to analyze. There- 
fore, for exposition purposes, they will be presented first. 
A detailed discussion for the justification of the choice of 
the boundary conditions will be presented elsewhere. 

In order to keep just a small number of parameters. 
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Eq. (2) neglects the force exerted by the magnetic field 
on the normal electrons. Therefore, the Hall voltage that 
we obtain is that of the superconducting electrons only. 
Since the Hall field is much smaller than the field in the 
direction of the current, the Hall voltage should be, to a 
good approximation, the superposition of contributions 
from both kinds of electrons. We will also assume that 
the density and mobility of the normal electrons are not 
significantly affected by superconductivity. For a treat- 
ment that considers the entire resistivity tensor see, e.g., 
Ref. 10. 

We integrate the TDGL equations by means of a finite 
difference method, using essentially the same program 
as in Ref. 8. In this method the sample is represented 
by a rectangular grid, consisting of x Ny cells with 
spacings a,j. — LlN,^. and ay = d/Ny. Discrete values 
of tp are defined at every vertex and values of A^^y are 
defined at every link in the respective x ov y direction. 
When periodic boundary conditions are used, they are 
imposed both on ip and on exp(«Aj,aj,). A standard ini- 
tial state was obtained by raising the applied field from 
to 0.5Hc2 and then keeping the field fixed until a stable 
(or metastablc) state was reached. The initial state for 
every run was then obtained by adding to ip at every ver- 
tex a complex random number with normal distribution, 
zero average and standard deviation 0.1. In a few cases, 
different histories were used. 

The current enters the algorithm through the effect it 
produces: it raises the value of the magnetic field at one 
interface and lowers it at the other. 

Most of the numerical studies in TDGL consider f] > 1, 
corresponding to cases in which TDGL can be derived 
from microscopic models; we have found that the eflFects 
reported here appear for rj significantly smaller than 1, 
as is the case for clean superconductors. 



III. RESULTS 
A. Periodic boundary conditions 

We report on samples of thickness d = 8^(0) « 7^(T). 
For these samples the standard initial state had a row 
of vortices at y = d/2 and the average distance between 
consecutive vortices was 2.5^(0) w 2.2^(T). 

Figure 1 shows the magnetization M of the film as a 
function of the current /, for several values of rj. The 
curves start at / = 4 x 10""^ rather than / = in order 
to chop off the influence of the random numbers added 
to the initial state. The general behavior is similar for 
all the curves (including additional values of r] not shown 
in the flgure). For every curve there is a small region 
< 7 < 7i where M is a smooth function of I. For 
/ > Ji, there arc points where the slope changes discon- 
tinously; at some of these points the slope changes sign 
and sometimes there is just a kink. (The curves them- 
selves have to be continuous if / changes at a finite rate.) 
The points where the slope is discontinuous coincide with 
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FIG. 1: Clusters of curves of the magnetization of the film 
versus the current, with as a parameter. Each curve in 

a cluster has an initial state slightly different than those of 
the others. The magnetization is the volume average of the 
induced magnetic field, in units of -H'c2(0), divided by 47r. 
The unit of current is c$o/[27r^(0)]^ per cm of length in the z- 
direction. For each cluster there is a different value of 77, which 
is marked next to it. Typically, each cluster contains four 
curves. For visibility, most clusters have been shifted in the 
vertical direction. Each curve starts at the current 4 x 10""^. 
At the right extreme of every cluster the film is in the normal 
state and the magnetization vanishes. The other parameters 
used in the calculations axe: = 80, Ny = 16, ax = ay = 
0.5, K = 2, T = 0.5 and the size of the noise is the same as in 
Rcfs. 7,8. The current increment between consecutive steps is 
AJ = 5 X 10^^ and the time increment varies from At = 0.008 
for 77 = 0.3 to At = 0.0015 for r; = 0.03. The applied field 
was kept fixed at 0.57/c2(0). In the inset, one of the curves 
for 77 = 0.03 is compared with the values obtained when the 
side of the cells in the calculation grid is decreased by a factor 
of 2 and the initial state is significantly different. 

the entrance or exit of vortices in and out of the film. 

The precise points where these discontinuities occur vary 
among different runs, but the general behavior is always 
the same. There is a current I2 > h where the magne- 
tization changes from diamagnetic to paramagnetic. If 
we disregard the rapid oscillations of the curves and con- 
sider only their smoothed trends, we observe that there is 
a current I3 > I2 where the slope decreases pronouncedly. 
(For T] = 0.03, I3 = 0.31.) Finally, there is a current Ici 
where the film becomes normal and the magnetization 
drops to zero. We also observe that there are regions, 
like a region that contains I = I3 for rj = 0.07, where os- 
cillations are practically absent and all the curves in the 
cluster coalesce. One might be tempted to suspect that 
in this regime vortices do not enter or leave the film, but 
closer examination shows that large regions of the sam- 
ple have become normal (typically, \tp\ < 10~^ in these 
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FIG. 2: M{I) for rj — 0.03. One curve is for increasing cur- 
rent, as in Fig. 1, and four curves are for decreasing current. 
The other parameters are as in Fig. 1. 



regions); under this condition, the concept of vorticity 
loses its significance. 

/c| is the Umit of metastability of the superconducting 
state, but at these high currents also the normal state 
is metastable. Figure 2 shows a cluster of M{T) curves 
for a case in which the current was decreased from above 
/cT to 0. We see that there is a limit of metastability 
Ici ^ 0.98 and the sample switches to the superconduct- 
ing state when this limit is approached. The curves M(/) 
are roughly the same when / is raised or lowered, except 
for a large hysteresis loop between I^,^ and I^i, and a 
small hysteresis loop between I3 and I31 ~ 0.24. All the 
curves in the cluster undergo a drop in the magnetiza- 
tion at the same current I^^. Hysteresis is a common 
phenomenon in superconductivity when the current is 
varied; a particularly well known case is that of an un- 
derdamped Josephson junction. 

In order to test the reproducibility of our results, we 
performed several runs with different computational pa- 
rameters. The solid line in the inset of Fig. 1 is one of 
the lines in the cluster for = 0.03, whereas the dots cor- 
respond to essentially the same physical parameters, but 
the computing grid was denser (A^^, = 159, A^^ = 32) and 
the rate of change of the current was increased by a fac- 
tor of 5 (At = 0.0003). In addition, the initial state was 
significantly different from the standard one; the average 
distance between consecutive vortices was 2.2^(0) rather 
than 2.5^(0). Due to the different initial state, the initial 
magnetization for this exceptional run is about half that 
of the runs in the cluster; however, for / > /a, it appears 
that the memory about the initial state has been lost. 

Figure 3 shows the average electric field E'y in the di- 
rection of the current, for a few values of rj. These results 
were obtained by averaging the parallel component of the 
electric field over all the lines in the grid in the direction 
of the current, and also over 2000 consecutive time steps. 
The results resemble those obtained for long channels in 
the absence of applied field. For 7 < 7i, E^^ is very small 
and for / > /c| (normal state) wc obtain = 2k^I /d, 
which is Ohm's law in the units we are using. In the 
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FIG. 3: Clusters of curves of the parallel component of the 
electric field versus the current, with 77 as a parameter. The 
electric field unit is c$o/[87r^K^C74^(0)], where a is the (nor- 
mal) conductivity of the sample. For visibility, the upper 
clusters have been shifted in the vertical direction. In this 
graph each hne starts at / = _E|| = 0. The other parameters 
used are the same as in Fig. 1. The inset shows the initial 
part of the considered curves. 



intermediate region, except for small oscillations, we ob- 
tain straight lines with the same slope as in the normal 
state. 

The inset in Fig. 3 is a close up for low currents. The 
reason that E\\ doesn't vanish completely for 7 < 7i is 
that our method of evaluation is not exactly stationary. 
For example, if we start from the point at 7 = 0.01 for 
T) = 0.07 and keep the current fixed, S|| decays with a 
time constant of 18 time units. This result can be un- 
derstood in terms of vortex motion: the current exerts 
a force on the vortices, which attempts to drive vortices 
into or out of the film. However, for 7 < 7i, the Bean- 
Livingston barrier^'^ prevents vortices from crossing the 
interface. Therefore, for constant currents the vortices 
attain equilibrium positions and stop moving; it is only 
the change in current that keeps the vortices in motion. 
Likewise, E\\ increases close to I\\ this happens because 
the configuration becomes unstable and vortices acceler- 
ate. 

Figure 4 presents the average of the y-component of 
the electric field, E±, for several values of rj. This re- 
sult was obtained by averaging this component over all 
the lines in the grid in the direction perpendicular to the 
current, and also over 2000 consecutive time steps. There 
are essentially two regions of currents for which E±_ does 
not vanish. There is a clearly distinguished feature at the 
right, present for 0.01 < rj < 0.07, which we call a "bub- 
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FIG. 4: Clusters of curves of the Hall component of the elec- 
tric field versus the current, with 77 as a parameter. Most 
clusters have been shifted in the vertical and/or horizontal 
direction, but they are all in the same scale. Each cluster 

starts at / = 4 X 10^'^ and at the right extreme E±_ = 0. The 
other parameters used are the same as in Fig. 1, except for 
A/, which was taken as 10"'^ for 77 = 0.01 and rj = 0.005, 
and as 2.5 x 10~^ for 77 = 0.15. E± usually vanishes, but for 
some regions of I E± ^ 0. In these regions E± is chaotic: 
for some runs is positive and for others is negative. For 
0.01 < r/ < 0.07 there is a region where E±{I) > for some 
lines and E±{I) < for others, but \E±{I)\ is the same for 
all the lines in the cluster, so that a "bubble" is formed. Note 
that the bubble for ij = 0.03 invades the neighboring clusters. 



ble"; there is also a less clear feature at the left, which 
we call a "jitter." In these regions E_\_ is chaotic: minute 
differences in the state of the system when it enters the 
region cause it to assume any of several very different 
functions E±{I). A major difference between these two 
regions is that for the bubble the system chooses among 
a small number of possibilities, which appear to bo sym- 
metric with respect to the line E±{I) = 0, whereas for 
the jitter the number of possibilities is large and has no 
obvious pattern. In most cases E± (I) is either positive or 
negative in the bubble region, but there arc a few cases 
(as a case shown for rj = 0.05) in which E±{T) = 0. 
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FIG. 5: Hall component of the electric field for different histo- 
ries. The solid line and the dots describe the same runs as in 
the inset of Fig. 1. In the absence of scattering, the dots look 
like a thick line. The dashed line describes a case in which 
the current was decreased. 



For 1] ~ 0.01 the bubble ends at /c|, but for 77 > 0.03 
the bubble ends at a current smaller than ■ In the case 
77 = 0.03, the bubble region coincides with the central 
part of the region H < I < Ic^ in which the magnetiza- 
tion curve looks lower and smoother than its continuation 
at both sides. As 77 increases, the bubble moves to the 
left, until it merges with the jitter. 

In order to elucidate which part of the film gives rise 
to E_\_, we have also evaluated its averages restricted to 
the sides of the cells that touch one of the interfaces. 
We found that E±^ in the bubble region is not influenced 
by the interface where the applied magnetic field is aug- 
mented by the induced magnetic field; E±_ in the jitter 
region is influenced by both interfaces. 

Since the Hall voltage we find might just be due to 
some instability of our numeric algorithm, we compare 
in Fig. 5 the values of E±_ obtained in different ways. 
Two curves are for currents increasing with time, but for 
grids and time steps of different sizes, and very different 
initial states (the same runs as in the inset of Fig. 1). 
Except for the amount of scattering and overshoot, the 
results coincide in the bubble region (although they have 
opposite signs in the jitter region). The third curve is for 
decreasing current. In spite of the large hysteresis found 
in Fig. 2, this line also coincides with those of increasing 
current in the region where the Hall voltage is present. 
It should be mentioned that for decreasing current E_\_ 
showed up in the bubble region only for about half of the 
runs. When it did show up, it appeared immediately with 
the switch to the superconducting state and remained 
present down to / = Izi- No visible change in E±_ is 
obtained if our boundary condition at y = and y = d 
is replaced by y ■ (— iV — A) ip = 0. 

We would like to gain some intuition concerning the 
reason for the existence of this chaotic Hall voltage. For 
this purpose, we have mapped the size of the order pa- 
rameter for several currents. Figure 6(a) is a contour 



FIG. 6: Contour plot of the size of the order parameter tp- W 
is larger in the lighter areas. The current flows in the vertical 
direction. The unit of length is ^(0). Note that the scale is 
different for each direction. All the parameters are those of 
the dotted line in the inset of Fig. 1. (a) / = 0.555, slightly 
before the bubble, (b) I = 1.195, almost at the end of the 
bubble, (c) / = 1.235, slightly above the bubble region. 



plot of IV-"! for a current slightly below the bubble region. 
This current is already large enough to turn into normal 
the entire area close to the interface where the magnetic 
field is large. The most interesting feature is that the 
superconducting region is not just a stripe parallel to 
the current direction, but it rather concentrates into a 
discrete set of superconducting islands; for the current 
in Figure 6(a), there are two such islands. Instead, for 
I = 1.235, slightly above the bubble region, we find that 
there are three well defined islands. In the bubble region 
itself, a process occurs in which two islands have to turn 
into three; during this process the islands assume an ir- 
regular shape, as in Fig. 6(b). This deformation imposes 
an overall motion of the superconducting part either in 
the direction of the current or against it, and this motion 
produces the Hall voltage. 

The size of the Hall voltage in the bubble region does 
not depend on the rate at which the current is swept. 
Increasing this rate by a factor of 40 just produces some 
extra overshooting, but near its maximum |i?i(/)| looks 
essentially the same as in Fig. 5. Moreover, if we stop 
increasing / and keep it at a fixed value in the bubble 
region, E±{I) remains constant in time, except for small 
oscillations that might be due to numeric inaccuracy. Ac- 
cording to our interpretation, this means that when I is 
in the bubble region the islands have to move, either in 
or against the direction of the current, even if I is kept 
unchanged. 

It is easy to monitor this motion by following the time 
dependence of E± . The effect of the motion is not ob- 
served if we take the average of E± along the entire range 
< a; < L, but becomes visible if we divide the sample 
into fringes, perpendicular to the current direction. Fig- 
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FIG. 7: Contour plot of E±{x,t). The unit of time is 
47rK^(T^(0)^/c^. In order to evaluate E±{x,t), the sample was 
divided into 8 fringes, perpendicular to the current, and E± 
was averaged on each fringe. For the case described here, 
r] = 0.03, / = 0.85, and the other parameters are as in Fig. 1. 
The regions x < 2.5 and x > 37.5 have been chopped off, due 
to the finite width of the fringes. This contour plot indicates 
that E±{x,t) behaves as a one-dimensional wave, with wave 
velocity 0.7cV47rH;VC(0). 



ure 7 is a contour plot for E±{x,t), which shows the 
motion of the maxima and minima of E±. Since these 
maxima and minima are expected to move together with 
the islands, their velocity should be equal to that of the 
islands. For the case described in Fig. 7, the islands move 
with a velocity of 0.7c^/47rK;^fT^(0). 

A quantity that is easier to interpret is the "number 
of superconducting electrons." By defining n^{x,t) = 
J\tjj{x,y,t)\'^dy and n^(y,t) = J \tp{x,y,t)\'^ dx and 
drawing contour plots of these quantities, we can visu- 
alize how the superconducting regions move. We have 
applied this procedure for the case tj = 0.03 while the 
current is kept constant at 7 = 0.4; for this current E± 
vanishes on the average. The upper panel in Fig. 8 shows 
that the superconducting regions do not move in the 
y-direction. The lower panel should be compared with 
Fig. 6(a): while in Fig. 6(a) there are only two islands, 
the distribution of exhibits four fringes. (Half a fringe 
was chopped off by the binning.) Examination of this 
lower panel shows that the islands switch positions with 
a period of about two time units, i.e., the centers of the 
islands become the valleys that separate between them, 
and vice versa. There is an additional phenomenon, with 
a period of about 8 units, in which the islands are blurred. 
In view of Fig. 8, Fig. 6(a) should be understood as an 
average over a period of time that is neither too long nor 
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FIG. 8: Distributions of the superconducting electrons in the 
X- and {/-directions, as functions of time, rj — 0.03, / = 0.4, 
and the other parameters are as in Fig. 1. The size of the 
binnings for x, t/ and t are 5, 1 and 1, respectively. 



a multiple of these periodic processes. 

From Fig. 8 we also learn that | should not be inter- 
preted as due to motion of the islands in the y-dircction. 
Rather, it acts like a Josephson field that causes oscil- 
latory motion of the superconducting regions in its own 
direction. When the volume average of E\ | is analyzed as 
a function of time, we find that is almost constant, with 
a small alternating part with period 8. 



B. Dirichlet condition 

In the previous case the islands formed a perfectly pe- 
riodic structure, due to the artificial requirement of peri- 
odic boundary conditions in the direction of the current. 
We may still anticipate that for any reasonable boundary 
condition at the electrodes there will be a finite number 
of islands and the creation of a new island will involve 
distortion. As a more realistic boundary condition at 
the electrodes, we considcircd Dirichlet conditions. The 
boundary condition for A was the assumption that the 
electric field in the electrodes is always in the x-direction, 
implying that Ay remains fixed in time. 

The results are complicated by the fact that now the 
electrodes pin the superconducting islands and impede 
their free motion. For L = 40^(0) the influence of pin- 
ning is so strong that we did not find island fragmenta- 
tion. We may expect that the influence of pinning will 
be weaker for longer samples, but, if the size of compu- 
tational cells is kept smaller than the coherence length. 
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FIG. 9: Magnetization (curve that looks like a mountain in 
the horizon) and Hall component of the electric field (curve 
that oscillates strongly about zero), for Dirichlet boundary 
conditions. Only one run was caxried out. The size of the 
sample is L = 80^(0), d = 16^(0). Other parameters: rj = 
0.03, K = 2,T = 0.5, aa, = ay = 0.5^(0), AI = 10"®, At = 
0.0015. 




(a) (b) (c) 




(d) (e) (f) 

FIG. 10: Contour plot of the size of the order parameter 
for the process described in Fig. 9. (a) / = 0.08. (b) I = 0.75. 
(c) 7 = 1. (d) 7 = 1.2. (e) 7 = 1.3. (f) 7 = 1.5. 
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larger L implies heavier numeric calculations. As a com- 
promise, we studied cases where L — 80^(0); we also 
doubled the thickness of the sample, although that was 
not really necessary. For these thicker samples the initial 
state had three rows of vortices. 

Figure 9 shows the magnetization and the average Hall 
field as functions of the current. Wc still obtain that E±_ 
vanishes for small and large currents and is present in 
an intermediate region. In this region there are two sub- 
regions where E±_ is significant, separated by a "quiet" 
subregion close to 7 ~ 1. For 1-05 </ < 1.45, E±_ has the 
same order of magnitude as in the "bubble" for periodic 
boundary conditions, but is now of the "jitter" type. 

Figure 10 shows the shape of the order parameter as 
the current is increased. For / = 0.08 there arc three 
rows of vortices in the sample, but around the row at 
the right il) is too small to make them visible. For this 
row the value of j-i/'l at the saddle points is of the order 
of 0.02. In a rough sense, there is a symmetry mirror at 
X = L/2 and, accordingly, we see that E±_ vanishes. For 
/ = 0.75 it would be useless to look at vortices; the entire 
area at the right has become normal and the notion of su- 
perconducting islands becomes more meaningful. There 
is no symmetry with respect to a; = L/2 and, accord- 
ingly, E± is large. For 7 = 1 a clear pattern of islands 
has developed. Two border islands arc pinned at the 
electrodes and four islands in the middle form a nearly 
periodic pattern. This pattern is symmetric with respect 
to X = L/2 and E± is small. For 7 = 1.5 we have again 
a regular symmetric pattern of islands and E± = 0. In 
contrast with the previous subsection, "counting" islands 
is not clear, since the border islands may be regarded as 
fractions that do not remain fixed. The cases 7 = 1.2 
and 7 = 1.3 illustrate the passage between the regular 
situations at 7 = 1 and 7 = 1.5. During this passage the 
regions where ip is significant move back and forth and 
assume irregular shapes. Accordingly, E±_ is large and 
changes sign. For large currents the total magnetic field 
is large at both superconducting-insulating boundaries 
(with opposite signs). Figure 10(f) shows that the super- 
conducting islands take refuge along a stripe where the 
magnetic field is small, and this enables them to survive 
up to large currents. 

Contrary to the previous subsection, the islands are 
now blocked by the electrodes and cannot move freely 
in or against the direction of the current; therefore, E±_ 
roughly vanishes on the average. However, islands might 
conceivably be formed near the middle of the sample and 
could migrate to the electrodes and disappear at them. 
If this is the case, we would expect E±{x) to differ from 
zero when averaged over fringes of limited width. Fig- 
ure 11 shows the averages of E± over fringes of width 
Ax = 5^(0). We see that E±{x) is roughly antisymmet- 
ric with respect to the line x = 40^(0) at the middle of 
the sample. We also see that E±{x) does not vanish in 
regions where the average E± does (7 < 0.1 and 7 > 1.5). 
The difference between these and the other regions is de- 
termined by whether the antisymmetry with respect to 




FIG. 11: The sample described in Fig. 9 was divided into 
16 fringes, perpendiculax to the current. The curves in the 
graph describe the average of E± in each of these fringes. 
The number next to some of the curves is the value of x 
at the middle of the fringe. The curves for fringes close to 
X = 40^(0) are hard to follow and have not been marked. 



the middle of the sample is perfect or not. In the fol- 
lowing subsection we will check whether the conjecture 
of island migration describes an appropriate scenario. 

In the regions where the average E± does not vanish 
and the islands are irregular, E±{x) is chaotic. Since this 
effect is very similar to the case that will be discussed in 
the following subsection, we do not provide a figure for 
it. 



C. Nonuniform sample 

It turns out that the results are more clear cut if the 
electrodes are "smeared" by gradually reducing super- 
conductivity in their proximity. 

A possibility for describing locally stronger or weaker 
superconducting materials is the replacement of Eq. (1) 
by 

dt^l^ = -- - A)' V + (1 - T) -1-6)^] +f , 

(3) 

where (5 is a function of position. If (5 > (respectively 
(5 < 0) in some place, superconductivity is stronger (re- 
spectively weaker) at that place; the case S = —1 de- 
scribes the situation in which the critical temperature 
has been reached. 

We considered a sample of thickness d = 8^(0) and 
length L = 120^(0). Its central segment 20 < a; < 100 
was uniform ((5 = 0), but close to the electrodes (0 < a; < 
20 and 100 < a; < 120) 5 varied linearly with x, reaching 
5 = —1 at the electrodes. In this way, the normal ma- 
terial at the electrodes was met when superconductivity 
had already disappeared. 

Figure 12 shows the results for the magnetization and 
for the average Hall field. Comparison of this figure with 
Fig. 9 shows that we have recovered some of the features 
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FIG. 12: Magnetization and average Hall component of the 
electric field for the sample considered in subsection IIIC. 
The averages of M and E_i_ were taken over the uniform cen- 
tral part (20 < a; < 100) only. Other parameters: d = 8^(0), 
rf = 0.03, n = 2,T = 0.5, = ay = 0.5^(0), A7 = 8 x lO"'', 
M = 0.0015. 
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FIG. 13: Hall field for difi'erent positions in the uniform seg- 
ment of a nonuniform sample. E±_ has been averaged over 
fringes of width 5^(0). Each curve is marked by the distance 
X between the middle of the fringe and one of the electrodes. 
For X > 60, can be obtained from E±{x) « — _Ex(120 — a;). 
For visibility, the results were grouped into pairs and the up- 
per pairs were shifted in the vertical direction. At / = 2, E± 
vanishes for all fringes. All parameters are the same as in 
Fig. 12. 



encountered for the case of periodic boundary conditions: 
there is an intermediate region where the Hall field van- 
ishes and there is only a limited region where this field is 
large. Also, there is a current where the over-all slope of 
the magnetization curve decreases significantly. On the 
other hand, there is no extended region where the Hall 
field preserves its sign. 

As in the previous subsection, the average Hall field 
is just a residual effect, due to its antisymmetry as a 
function of position. In Fig. 13 wc present E±_ for sev- 
eral values of x within the segment where the sample is 
uniform. Since antisymmetry is quite well obeyed, these 
results are shown for half of the sample only. 

In order to check whether E±{x) is chaotic, we re- 



FIG. 14: Hall field averaged over fringes of width 5^(0) for 
miscellaneous values of x and for a range of currents that 
contains the region where the average over x is large. All 
parameters are the same as in Fig. 12, but there are clusters 
of curves starting at four different initial states. For visibility, 
the cluster for x = 77.5 has been shifted by 0.05 to the right 
and 0.0012 upwards; the cluster for x = 97.5 has been shifted 
by 0.15 to the right and also 0.0012 upwards. 



120 




FIG. 15: Contour plots for the size of the order parameter for 
I = 0.88 (left) and I = 1.76 (right). The other parameters 
are as in Fig. 12. 



peated the calculation of Fig. 13 for several initial states. 

Some of the clusters of curves obtained arc shown in 
Fig. 14. We see that in the regions where the average 
value of E± is small in Fig. 12 all curves practically co- 
alesce, whereas in the region 1.25 ^ x < 1.5 there is 
chaotic behavior. 

As in the previous cases, we would like to relate the 
behavior of the Hall field to the shape of \tp\. Figure 15 
shows these shapes for currents where the average Hall 
field vanishes, before and after the region where this 
average is large. As might have been expected, these 
shapes are symmetric with respect to the transformation 
X 120 — X. A striking feature is that for J = 1.76 su- 
perconductivity concentrates into two sharply bounded 
islands at the borders between the segment where the 
material is uniform and those where superconductivity 
gradually decreases. These islands may be regarded as a 
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FIG. 16: Contour plots for the average density of supercon- 
ducting electrons along fringes perpendicular to the current, 
for the sample considered in subsection IIIC. For distances 
smaller than 15^(0) from an electrode, |-i/)(a;,t)p is very small 
and these regions are not shown. Left: I = 0.88. The size of 
the binning along the t-axis is 1, and along the x-axis is 5. 
Right: / = 1.32. The size of the binning along the t-axis is 2, 
and along the x-axis is 5. 

generalization of surface superconductivity and are prob- 
ably the reason for the long tail of the magnetization in 
Fig. 12. 

We now attempt to find a relation between E± and 
island motion. As in the lower panel of Fig. 8, the contour 

plots in Fig. 16 describe the density of superconducting 
pairs, \t{j{x,t)\'^ , averaged over fringes of width Ax = 
5^(0). The plot at the left is for a current such that 
E± = in Fig. 12 and at the right for a current such that 
E± is large. The plot at the left looks like a standing wave 
superimposed on a nonuniform background. This means 
that there is no net motion of the islands, and the fact 
that E±(x, t) does not vanish for every value of x should 
be attributed to a nonlinear effect. The plot at the right 
might be described as a case of "breathing." The most 
notable feature of this plot is that \ip{x, t)\ decreases and 
increases periodically for .x ^ 65. In addition, wc see that 
the island-pattern moves asymmetrically: there is no net 
motion for a; < 50, whereas there is motion away from 
the center for x > 70. This confirms the conjecture that 
the existence of an average Hall field is related to island 
motion. 



IV. DISCUSSION 

We have found that for appropriate parameters a 
nearly periodic pattern of superconducting islands can 
form in a superconducting sample in the Hall configura- 
tion. There are situations for which a pattern becomes 
unstable and the islands become irregular, until a new 
regular pattern is attained. While the islands are ar- 
ranged irregularly, they move in (or against) the direction 
of the current and the contribution of the superconduct- 
ing electrons to the average Hall voltage does not vanish. 



Even when the islands do not move in a preferred direc- 
tion, the size of the order parameter oscillates in time, so 
that local maxima of the order parameter may periodi- 
cally become local minima. 

This scenario shows up very clearly when periodic 
boundary conditions are imposed at the electrodes. For 
other boundary conditions, these phenomena are still 
qualitatively present, but become more difficult to an- 
alyze. The difference between realistic and periodic 
boundary conditions is that the latter don't pin the su- 
perconducing islands. 

We might mimic a situation in which islands are not 
pinned by taking a sample which is infinitely long in the 
direction of the current flow; however, we may expect 
that no Hall voltage would appear in this limit, since 
the islands structure could shrink gradually and no dis- 
tortion would be required. It seems therefore that some 
amount of pinning is desirable. This is consistent with 
the theoretical results obtained in Ref. 10: the strongest 
oscillations of the Hall voltage were obtained for an inter- 
mediate amount of defects. In the calculations of Ref. 10 
the defects were assumed to be located in a periodic ar- 
ray: this arrangement should be expected to favor the 
formation of periodic arrays of islands. Artificially cre- 
ated periodic arrays of defects in superconductors have 
been available for a long time.'^'* We also found that grad- 
ual decrease of the superconducting strength of the su- 
perconducting material near the electrodes may help to 
create conditions similar to the periodic case. 

One might argue that there is no qualitative difference 
between the formation of a new island and the entrance 
of a vortex: when a vortex enters the sample, the oth- 
ers have to move away from it and this motion gives rise 
to Hall voltage. There is, however, a quantitative differ- 
ence: vortices feel each other over distances of the order 
of the coherence length, whereas islands seem to feel each 
other over distances that are larger by at least an order 
of magnitude. 

We do not have a tested hypothesis concerning the 
physical parameters that control the size of the islands. 
Conceivably, the length of the islands could be propor- 
tional to their width and this could be controlled by the 
width of the region where the magnetic field is sufficiently 
small. On the other hand, since the effects described in 
this article have only been found for small values of r], 
we might suspect that the size of an island is related to 
the penetration length of the electric field, ^{T)/y^} 

Our results are indirectly supported by the experi- 
ments in Ref. 10. In that reference, which used a par- 
ticularly clean material (so that r] was indeed expected 
to be small), the variable quantity was the applied field, 
whereas the current was fixed. However, the Hall voltage 
exhibits the same qualitative features as in our case: it is 
present in limited regions only, is absent close to the crit- 
ical current, exhibits oscillations and its sign is reversed 
for different experiments that were apparently identical. 

The fact that we do obtain a Hall voltage using Eq. (1) 
challenges the accepted contention that a transverse 
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voltage cannot arise from a purely dissipative TDGL 
model. ^'^^ Note however that the voltage we have found 
is statistically symmetric about zero, so that our result 
may be regarded as a sort of symmetry breaking. 

It should be possible to use some of the techniques that 
are employed in visualization of vortices (such as STM) 
for the visualization of the islands predicted here. 
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